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Abstract. This paper deals with the early results of a new model of pedestrian flow, conceived 
within a measure-theoretical framework. The modeling approach consists in a discrete-time 
Eulerian macroscopic representation of the system via a family of measures which, pushed 
forward by some motion mappings, provide an estimate of the space occupancy by pedestrians 
at successive time steps. From the modeling point of view, this setting is particularly suitable 
to treat nonlocal interactions among pedestrians, obstacles, and wall boundary conditions. In 
addition, analysis and numerical approximation of the resulting mathematical structures, which 
is the main target of this work, follow more easily and straightforwardly than in case of standard 
hyperbolic conservation laws, also used in the specialized literature by some Authors to address 
analogous problems. 



1. Introduction 

In the last decades, the modeling of crowd motion and pedestrian flow has drawn the attention 
of applied mathematicians, because of an increasing number of applications, in engineering and 
social sciences, dealing with this or similar complex systems. First studies on pedestrian behavior, 
dating back to the sixties, were aimed at providing guidelines for the design of walkway infras- 
tructures in urban areas, and yet nowadays crowd-structure interaction is the object of scientific 
investigations (see e.g., Venuti et al. |17j and the main references therein). More recently, man- 
agement and optimization of pedestrian fluxes in large and crowded environments, like airports, 
stations, shopping centers, stadiums, have motivated new efforts toward more accurate experi- 
mental investigations (see e.g., Helbing et al. [HIHD]) and a more targeted mathematical modeling 
of human walking attitudes in bounded areas, with special emphasis on the concepts of "walking 
program" of pedestrians and of pedestrian-pedestrian, pedestrian-obstacle interactions. 

The existing literature on mathematical modeling of human crowds is especially concerned with 
models at the microscopic scale, in which pedestrians are described individually in their motion 
by a set of ordinary differential equations. Problems are usually set in two-dimensional domains 
delimiting the walking area under consideration, and the presence of obstacles within the domain, 
as well as of targets that pedestrians aim at, is possibly taken into account. The basic modeling 
framework relies on classical Newtonian laws of point mechanics: Pedestrians are assimilated to 
rigid disks of fixed, possibly different, radii, with a velocity, or alternatively an acceleration, that 
takes into account both the desired direction of motion and a superimposed perturbation due to 
the presence of other pedestrians in the surroundings (Helbing and Johansson [7 , Maury and Venel 
[15)). In practice, pedestrians are assumed to move mainly according to a desired velocity, i.e., the 
velocity they would have in the absence of other people, so as to get a desired target and avoid at 
the same time intermediate obstacles. The desired velocity is however modified by the necessity to 
fulfill some local geometrical constraints of maximal congestion. More complicated strategies to 
determine the direction of motion may involve the minimization of some 'walking cost', that each 
pedestrian evaluates on the basis of the current configuration of the system in her/his neighborhood 
and of her/his prediction about the motion of other pedestrians (Hoogendoorn and Bovy |12j). 

A different approach to the problem, however not yet as much developed in the specialized 
literature, consists in using partial differential equations at the macroscopic scale, that is in de- 
scribing the evolution in time and space of pedestrian density rather than following each subject 
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individually. The starting point is the analogy with classical fluid dynamics: An Eulerian point 
of view is adopted, appealing to the conservation of mass of pedestrians supplemented by either 
suitable closure relations linking the velocity of the latter to their density (see e.g., Hughes |13j ) 
or an analogous balance law for the momentum (see e.g., Bellomo and Dogbe [2]). Again, typical 
guidelines in devising this kind of models are the concepts of preferred direction of motion and 
discomfort at high densities. In the framework of scalar conservation laws, a macroscopic one- 
dimensional model has been proposed by Colombo and Rosini [5] , resorting to some common ideas 
to vehicular traffic modeling (e.g., fundamental diagrams), with the specific aim of describing the 
transition from normal to panic conditions. 

Microscopic models allow an accurate description of the behavior of each single agent of the 
system, but they usually require a large number of ordinary differential equations to be dealt with 
simultaneously, as many as the number of pedestrians considered in the model. This often makes 
it difficult, for practical purposes, to recover macroscopic information on the system from the 
knowledge of the microscopic state, hence to investigate the relevant global features possibly also 
in connection with control and optimization problems. On the other hand, macroscopic models 
existing in the literature heavily rely on the fluid dynamics analogy, hence they use mathematical 
tools proper of hyperbolic conservation laws, which however might not be the most appropriate 
ones to address the behavior of pedestrians. For instance, hyperbolic equations require the impo- 
sition of boundary conditions in a technically tricky way, which may not correspond to the actual 
modeling needs. Moreover, crowd motion is definitely not one-dimensional: A two-dimensional 
description is inevitable in order to catch the fundamental aspects of the system, but it is known 
that multidimensional hyperbolic equations generate additional analytical and numerical difficul- 
ties, leading to an increase of technical complexity in handling the final product. 

Bearing all above in mind, we propose to adopt a different macroscopic point of view, based 
on a measure-theoretical framework which has recently been introduced by Canuto et al. @] for 
coordination problems (rendez-vous) of multiagent systems. 

Given a two-dimensional spatial domain fl C K 2 , possibly containing obstacles understood as 
internal boundaries, the basic idea is to describe the space occupancy by pedestrians at time n 
via a measure /x n , that is, roughly speaking, a mapping that to each subset E of f2 associates 
a real nonnegative number fj. n (E) representing an estimate, in macroscopic averaged terms, of 
the amount of people contained in E. The whole mass of pedestrians is then obtained as /U„(fi), 
and its conservation in time is possibly achieved by requiring /^ n (f2) = (j.q(£1) for all n > 0. If 
v n — v n(%) is the velocity field at time n of pedestrians located at a point x € f2, the dynamics of 
the system is described by a family of mappings {"f n }n>o such that J n { x ) — x = v n (x)At, which 
entails the following evolution for the measure fj, n (push forward) : 

Vn+l(E)=fl n ( 1 n 1 (E)), VECtl. 

This is nothing but a formal mathematical statement of the simple idea that the amount of people 
contained in a spatial region E at time n > is related to the analogous amount at the initial time 
n = along the trajectories of the motion of pedestrians themselves. In this context, we recall 
that a measure theoretical approach for mass transportation problems has been proposed also by 
Buttazzo et al. [3], focusing in particular on concentration and congestion effects. Specifically, 
the Authors claim that the latter may be applied to the modeling of crowd dynamics in panic 
situations. 

The construction of the velocity v n is the main modeling task. The framework we use easily 
allows to duly incorporate the concept of desired velocity, as well as to model interactions among 
pedestrians, possibly accounting for averaged non-local effects due to the surrounding crowding. 
By naturally tracing, for fixed x G Q, the displacements of pedestrians, the mapping j n (x) de- 
scribes the system under an essentially Lagrangian point of view, which is, in a sense, a more 
natural way to look at the motion of the agents. On the other hand, the push forward of the 
measure fj, n finally refers to an Eulerian handling of the system, and takes thus advantage of a 
fixed geometry without the need for resorting to the concept of evolving reference configuration. 
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The paper is organized into three more sections that follow this Introduction. In Sect. [2] the 
modeling framework by time-evolving measures is introduced for generic d-dimensional continu- 
ous systems. After a first parallelism with the classical mass conservation equation of continuum 
mechanics, the theory focuses on discrete time push forward, which is shown to represent an ex- 
plicit time discretization of a conservation law for a family of continuous-time-evolving measures. 
Well posedness and spatial approximation of the discrete time model are addressed. In addition, 
a computational scheme for its numerical treatment is derived, with a related convergence and 
error analysis. In Sect. [3] the above theoretical framework is specifically applied to pedestrian 
flow. From the modeling point of view, this setting proves to be useful in extending the idea of 
macroscopic mass to cases, like in fact human crowds, in which mass and mass density may not be 
conceptually defined in the straightforward sense of continuum mechanics. It may be questioned 
that this problem is actually common to many other classical and nonclassical systems, ranging 
from gas dynamics to vehicular traffic, for which, despite the essential granular nature of the mat- 
ter under consideration, the continuum approach has been successfully used. Nevertheless, in the 
case of pedestrian flow the measure theoretical approach allows to deal with nonlocal interactions 
among pedestrians, obstacles and more generally wall boundary conditions, and numerical ap- 
proximations more effectively than standard hyperbolic conservation laws. Finally, Sect. [4] draws 
some conclusions and briefly sketches research perspectives. 



2. Time-evolving measures theoretical framework 

Many real world systems can be mathematically described at the macroscopic scale, usually 
invoking conservation or balance laws of some gross quantities. As a matter of fact, in modern 
applications dealing with life and social sciences, e.g., biological systems, vehicular traffic, pedes- 
trian flow, often the only conservation principle which is reasonable to resort to is the conservation 
of mass, because the evolution of the system involves complex behavioral aspects which need not 
preserve other classical mechanical quantities, like linear momentum or energy. 

In the physical space K rf (d = 1, 2, 3 for applications), the mass can be viewed as a Radon 
positive measure p : B(M. d ) — > R + , i.e., a mapping from the Borel c-algebra B(M. d ) to the real 
nonnegative line, whose countable additivity furnishes the mathematical counterpart of the prin- 
ciple of additivity of the mass. The Lebesgue measure C d on W l formalizes instead the concept of 
volume. If one assumes the continuum hypothesis, stating that the mass p is absolutely continuous 
with respect to C d (written p <C C d ), namely that every body with zero volume has also zero 
mass, then Radon- Nikodym theorem implies the existence of a function p g Lj oc (M. d ), p > a.e. 
in M. d , called the moss density, such that dp, = p dx. 

In continuum mechanics, for a mass density pt evolving in timcj^jthe principle of conservation 
of the mass is stated in Lagrangian form as 



d 

dt 

r t (B) 



p t (x)dx = 0, V£e6(R d ), (2.1) 



where T t is the motion mapping, which describes the motion of the points of M. d under the action 
of a certain velocity field Vt as 

f t (a:) = V t (T t (x)), 

supplemented by the further condition Tq(x) = x. In practice, T t (x) is the position occupied at 



time t by the point initially located at x. Using Reynolds theorem, Eq. (2.1) can be formally 



converted in Eulcrian form and rewritten as the well-known conservation law for the density pt- 

^ + V-(p t F t ) = 0. (2.2) 



Throughout this section, the subscript t does not denote the partial derivative with respect to time. We use 
the notation pt(x), instead of the more classical p(t, x) (and similarly for other functions), to emphasize the idea 
that the density should be regarded as a function of x £ R d parameterized by t > 0. 
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Equivalence between Eqs. (2.1| and (2.2), however, holds for smooth p t and V t only, as passing 



from Eq. (2.1 1 to Eq. (2.2) requires some manipulations which are not valid if the involved 



functions are not smooth. 

As a natural generalization of Eq. (2.2) to the case in which a mass density need not be 
present, one can assume that a family of time-evolving measures fa '■ B(M. d ) — > K + , t > 0, is given, 
satisfying 

BlH ' « - ^ " (2.3) 



- + V • (ptV t ) = 0. 



dt 

This partial differential equation has to be meant in the sense of measures as follows: For every 
infinitely diffcrentiable test function r\ with compact support in K d , i.e., r\ £ Cfi°(M. d ), 

d 



dt 



r)(x)dfa(x) = / Vt](x) ■ Vt(x) dfa(x). 



(2.4) 



In particular, a family of measures {fa}t>o is said to be a solution to Eq. (2.3 1 if for all rj g (M. d ) 
the mapping 

t !->• / Tj(x) dfa(x) 



is absolutely continuous and satisfies Eq. (2.4). Notice that a basic requirement for the right-hand 



side of Eq. Q to be well defined is V t G (P(M d , fa)) d for all t > 0, i.e., 

\V t (x)\dfa(x) <+oo, V£>0. 



Let us denote by suppfa the support of the measure fa, and let us assume that an open set 



£1 C R d exists such that supp/^ CC for all t > 0. We can then rewrite Eq. (2.4) by integrating 
on fl only: 

j t J v{x)dn t (x) = J W V (x) ■ V t (x)diH(x), V77 e C °°(O). (2.5) 

SI SI 

In addition, there exist an open set U C f2 and a test function r\u g C q°(0 ) such that supp/x t C U 
for each t > 0, ry^j = 1 on supp^ t , and 77(7 = on f7 \ U, whence Eq. (2.5) with the choice rj = r\jj 



(2.6) 



gives 4f fa ( SU PP Mi) = 0- However, /x t (supp /x t ) = jti t (0) for all t > 0, thus finally 

which states that no matter is flowing through the boundary <9f2 at any time: The measure of Q 
is not varying, hence the whole mass remains concentrated within it. Condition supp fa CC f2 for 
all t > compares therefore with a classical homogeneous Neumann boundary condition for the 



PDE (2.2) 



2.1. Push forward and time discretization. When dealing with a discrete time evolution of 
the system, one can fix a time step At > and trace the mass by a sequence of positive measures 
{fai}n>o, fai ■ B(R d ) — > R + each n, via the following recurrence relation (push forward): 

fai+l = 7n#A*n, (2.7) 

where 7„ : R d -> R d is the one-step motion mapping (briefly termed motion mapping in the sequel 
for simplicity). More specifically, 

7 n (z) = x + v n (x)At, (2.8) 
with v n : R d — > R d a velocity field, so that "f n (x) is the position at the (n + l)-th time step of the 
point which at the n-th time step is located in x. 



The push forward (2.7) has to be understood formally as 

» n+1 (E)=fa l ( 1 - 1 (E)), V£e£(K d ), (2.9) 

which shows that fai+l is unaffected by the values that j n possibly takes outside supp/i„. On the 
other hand, supp// n _|_i = 7„(supp /i„), hence if there exists a set f2 C R d such that 7 n (f2) C for 
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is prescribed with supp //q C fi, then 



all n, and furthermore an initial measure /i : B(R d ) — > 
supp (J, n C il for all n > 0, whence we get, analogously to Eq. (2.6), the conservation of the mass 
of O: 

£tn(n) = ^oWi Vn>0. 



The definition of push forward given by Eq. (2.9) is equivalent to 

r)(x) dfx n+ i(x) = I rj{ln{x)) dn n (x) 



(2.10) 



for every bounded and Borel function r\ : M. d — > M. In particular, if 77 g C§°(M. d ) we can expand 
?7(7nO)) = V(x + v n (x)At) = r)(x) + AtVrj{x) ■ v n (x) + 0(At 2 ) 



and then plug this expression into Eq. (2.10) to get 



1 

At 



n{x) d/i„+i - / rj(x) d/i 



Vr)(x) ■ v n (x) dfx n + O(At), 



which, at least for /j„-uniformly bounded v n , can be viewed as an explicit time discretization of 



Eq. (2.4) up to identifying v n (x) = V n At(x). Again, if supp /i rl C fl for all n > then Eq. (2.10) 
can be rewritten by integrating on f2 only: 

T](x) dfi n+ i (x) = J ri(-f n (x))dfi n (x), (2.11) 
and the same Taylor expansion performed above provides now an explicit time discretization of 



Eq. (2.51 for test functions 77 <E C§°(0). 



Proposition 1. The push forward ( |2.7[ ) of the measure fj, n is the direct time discretization of the 
conservation law 

^u t (T t (E))=0, VEeB(R d ). (2.12) 
at 

Proof. To see this, discretize the motion mapping T t by the family of one-step motion mappings 
{7n}n>o as r„ + i = 7„ o • • • o 70, then approximate the time derivative in Eq. (2.12) above as 

(T n+1 (E)) - fi n {r n {E)) 
At 

Setting E :— T n+ i(E), and consequently r„(£') = 7 r] ; 1 (£ , ), yields 

fi n +i{E) = Unij^iE)), 



= 0. 



whence we recover formally the relation /i n +i = 7n#/-i n (cf. also Eq. (2.9)) or, equivalently, Eqs. 
plol), pTTTl). □ 



Remark. The Lagrangian mass conservation law (2.1) is a particular case of Eq. (2.12) for the 



measure d/j, t = p t dx. It is worth noticing that, unlike classical procedures to derive the pointwise 



Eulerian mass conservation equation (2.2), or analogously one of its weak forms (2.4), (2.5), from 



Eq. ( 2.1 ), the time discretization allows very few regularity of the involved fields in order to attain 



to Eqs. (2.10), (2.11) fromEq. (2.12) 



In order to deal with problems posed in a fixed bounded domain fi C R d , from now on we con- 
centrate on motion mappings j n : Q — > Q such that j n (£1) C f2. This way, given an initial measure 



Ho supported in fi, all measures /j, n deduced from the recurrence relation (2.7) are supported in 



f2 as well, hence we can simply think of them as defined on the Borel er-algebra B(Q) and refer 
directly to Eq. (2.11 1 whenever necessary. Notice that the measure of O is conserved in this case, 



indeed p n+ i(Jl) = /j„(7 ri ^ 1 (f2)) = fi n (Q) and the claim easily follows by induction. In particular, 
if /io(^) < +00 then fl will have a finite measure for all successive times n > 0. 
We begin by stating some basic properties of the motion mappings. 



Assumption 1. For each neN, we assume that the motion mapping 7„ 
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(1) is Borel, i.e., 7 ~ 1 (_E) <= B{fl) for all E £ B{tt); 

(2) satisfies 

^(^(EjjKCC^E), V£eB(fi), (2.13) 
for a certain constant C > 0. 



Let us briefly comment on property Q of Assumption [I] Given a Borel set BCfl, Eq. (2.131 



allows to control the Lebesgue measure of the inverse image of E by the Lebesgue measure of E 
itself. From a different point of view, it requires that 7„ does not map Lebesgue non-negligible 
subsets of SI into Lebesgue negligible sets. Indeed, if A C SI is such that C d (A) > 0, it is 



impossible to have C d {^ n {A)) = 0, for Eq. ( |2~13| implies C d (-f n (A)) > C- x L d {A) > (take 
formally E = j n (A)). 

Property ^ is for instance satisfied if the following holds true: 

|z-y|<A|7 n (a;)-7 n (y)|, V x, y E (2.14) 



for a certain constant A > 0. In this case, a suitable constant is C = A d , indeed Eq. (2.14) entails 



on the one hand £ d ( 7jl 1 (E)) < A d £ d ('y n ("f n 1 (E))), while on the other hand it suffices to observe 



that jnhnHE)) C E, thus £ d ( 7 „ (7- 1 (E))) < C d (E). We claim that the relation ( [2~14| holds if, 
for instance, the velocity v n is Lipschitz continuous on Q with Lipschitz constant < L < At : 

\v n (x) - v n (y)\ < L\x - y\, Vx,y€ti. 

In fact, in such a case we have 

\ln{x) - 7n(y)\ =\{x-y) + M{v n (x) - v n (y))\ 

> \x - y\ - At\v n (x) - v n {y)\ 

> (l-LAt)\x-y\, 



and if L fulfills the previous constraints we can take A = (1 — LAt) > 1. 
Another possible motion mapping complying with property ([2| is 

M 

7 „ (x) = x + At ^2 a iXEi (x), 

i=l 

where {E^fl-^ is a pairwise disjoint partition of SI, i.e., Int_E j; n IntEj = each i ^ j and 
Lif£ 1 Ei = SI, XEi is the characteristic function of Ei, and a™ € R d are constant, so that the 
velocity v n turns out to be the piecewise constant function v n (x) = J2i=i a ™XBi( x )- Such a 7„ is 
a measurable piecewise translation on S7 with the following property: For any E 6 S(S1), it results 

M M 

^(in'iE)) = Y J C d ( 1 - 1 {E) n Ek) ^£ d (J?n 7n (£,:)) < MC d {E), 



due to the invariance of Lebesgue measure under translations. A possible constant for Eq. (2.131 
is thus C = M, the number of elements of the partition, although it may not be the optimal (i.e., 
the smallest) one. For instance, in case of constant velocity v n (x) = a E R d in SI, which gives 
a linear transport "f n (x) = x + aAt\n{x) , the sets j n (Ei) are pairwise disjoint and the last sum 



in the computation above equals C d (E n 7„(f2)) < C d (E), so that Eq. (2.13) holds now more 
precisely with C = 1 regardless of M. 

If the initial measure no is absolutely continuous with respect to C d , a natural question is 
whether the same is true for all other measures fj, n recursively generated by the push forward. 
Property ^ of the motion mappings turns out to be designed precisely for this purpose, indeed 
we have: 

Theorem 2. If fi < C d then fj, n < C d for all n > 0. 

Proof. We proceed inductively on n. Assume fj. n <C C d for a certain n and consider E g £>(S1) 
such that C d (E) = 0. Then C d ( 1 - l {E)) < CC d (E) = 0, whence 

f , n+1 (E)= f i n ( 1 - 1 (E))=0, 

i.e., n n+ i <C C d . Owing to /i <C C d , we finally get by induction the thesis. □ 
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When a density p n £ ^(fl), p n > a.e. in Q, exists for the measures p n , Eq. (2.11) rewrites 

rj(x)p n+1 (x) dx = J r){~f n (x))p n {x) i 



) dx 



(2.15) 



for all bounded and Borel functions 77 : CI — > R. The new unknowns of the problem are now the 
p n 's, for which we can state the following properties. 

Theorem 3. Let po € L (f2) 6e given, po > a.e. in fi. T/ien i/iere exists a unique sequence 
{pn}n>o C L 1 (fi), p„ > a.e. in CI, solving Eq. (2.15) with po as initial datum, and moreover 

||Pn||i = ||po||i, Vn>0. 
If, in addition, p Q € L 1 (f2) n i/ien afeo p„ G L^fi) n L°°(f7) wrf/i 

llpnlloo < C n || Po ||oo, Vn>0. 



Proof. Existence of a sequence {p n } n >o C L (CI), p n > a.e. in CI, solving Eq. (2.15) is implied 
by Theorem [2] if one understands po as the density of po with respect to C d . 

(i) To obtain uniqueness, assume first that p n +i, p' n +i S L (CI) are such that 
n(x)p n+ i(x) dx = [ r](x)p' n+1 (x)dx = [ n(~/ n (x))p n (x) dx, 

Q 

that is 

??(a;)(p n+ i(x) - da; = 

for all bounded and Borel functions 77 : 17 — >• K. Confining the attention to the continuous 
77 compactly supported in CI, i.e., 77 g Cq(CI), we see that, owing to Riesz representation 
theorem, this relation defines the null functional in the dual space (Co(O))'. Therefore 
(Pn+i ~ Pn+i)dx must be the null measure, which implies p n+ i(x) = p' n+1 (x) for a.e. 
x G Cl, whence the ^-uniqueness of the density at the (n + l)-th time step. Proceeding 
now inductively from po, we get the uniqueness in L 1 (fi) of the sequence of densities p n . 

(ii) When estimating the i^-norm of the p n 's we can take advantage of their nonnegativity to 
discover, from Eq. (2.15) with i] = \n- 



||Pn+l||l = / Pn{x) dx = \\p n \\i. 

h 

Thus, by induction, ||p n ||i = ||po||i each n > as desired, 
(iii) In order to prove the last statement of the theorem, we consider a pairwise disjoint partition 
{Ei}f^ 1 of Cl and construct the following step function: 



M 



s ( x ) = ^2 a iXEi(x), Oti 



t = l 



C d {Ei) 



Pn{x) dx 



for every Ei such that C d (E t ) ^ 0. Notice that 



M M „ „ 

s||i = = ^ / p n (x)dx= I p n (x) dx = \\p n \\] 



and moreover CKj > 0, therefore s € L 1 (J7), s > a.e. in Cl. If C d (Ei) = for some i, 
then the corresponding coefficient on can be arbitrarily defined without affecting s as an 
element of 
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(iv) Set 



:= max C d (Ei) 



M 



We claim that, for 



> + , s converges pointwise almost everywhere in SI to p n +i- For 
this we observe, first of all, that, taking r\ = XE t in Eq. (2.15), the coefficients on can be 
rewritten as 



07 



Pn+l{x)dx - 



Ei 



a 



i{x) 



and 



Then we fix x G SI and notice that there exists a unique index 1 < i = i(x) < M such that 
x G -E^r,,:)) because the partition is pairwise disjoint. Thus s(x 

1 



\Pn+l{x) ~ S{X)\ 



Pn+l{x) 



C d {E l{x) 



< 



1 



C d {E i{x) ) 



Pn+i(y)dy 

Pn+l{x) ~ Pn+l{y)\ dy 



As m — > + , Ei( x } shrinks to x and Lebesgue differentiation theorem implies that the 
right-hand side of the above inequality tends to zero, therefore 

lim \p n+ i(x) — s(x)\ = for a.e. x G O, 

m-»0+ 

which gives the desired pointwise convergence, 
(v) Assume p„ G L X (Q) n then 



, £ d ( 7 - 1 (^)) , l .I .^1 
on < — „ Jr „, — HPnlloo < C\\p r , 



v< = 1, 



M, 



hence s G L°°(S1) with s(x) < CUpnHoo for all x G ft. Consequently, in view of the 
pointwise convergence, also p n +i(x) < C\\p n \\ oo for a.e. x G £1, which says p n +i G L°°(f2) 
with 

||Pn+l||oo < C||Pnl|oo- 

Proceeding inductively from po yields finally p n G L°°(f2) with the desired estimate on 
the L°°-norm. □ 

2.2. Spatial approximation. We turn now our attention to the construction of a spatial approx- 
imation of the push forward (2.151, which will result in a computational scheme for the numerical 
treatment of the problem 



'Find {p n }n>o C L 1 (f2) n L°°(d), p n > a.e. in SI, for all n > 0, such that 



v( x ) Pn+i(x) dx = I rj{"f n {x)) p n (x) dx, V?7 : O — > R bounded and Borel 



(2.16) 



for a prescribed initial datum p € L 1 ^) n L°°(f2), po > a.e. in f2. Theorem [3] gives the well 
posedness of this problem, i.e., existence and uniqueness of the solution, with suitable estimates 
on the L 1 and L°°-norms in terms of the corresponding norms of the initial datum. 

To be definite, we consider as domain SI a (hyper)cube in M. d , that we partition by a family of 
nested pairwise disjoint grids {-E?} 7= \ made in turn of (hyper)cubes with edges of constant length 



h > 0, so that C d = h d for each j: 



Mh 



\jE$ = Sl, Int^Jnlnt^) =0 Vi^j. 

The nesting of the grids implies that, given < h" < h' , for all < j < Mh" there exists 
< i < M h i such that Ej C E^ . Essentially the same construction described below may be in 
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principle repeated for domains and/or grid elements of different shapes, up to possible technicalities 
in the practical implementation of the resulting numerical scheme. 



piecewise constant functions over the grid {Ej}^^: 



One of the most natural ways to approximate Problem (|2.16|) is to introduce the following 

<h\M 
3 Jj = 



Ph{x) = Y d PlhX E ${x)> n>0, 
3=1 



such that P£(x) = p] h for all x € E%, and to replace p n , p n+1 in Eq. ( pl5| by Pg, P™ +1 
respectively. Notice that each P£ is integrable over f2, with 

M h 
3=1 

In addition, one can in general expect that the true function j n , which operates the push 
forward from the time step n to the time step n + 1, is not known exactly: Most frequently, an 
approximation g% : ft — > 51 of it is available, resulting from some discretization of the velocity v n . 
In the sequel, we will consider specifically 

gl(x)=x + ul(x)M, 

where u\ is piecewise constant over the grid {JSH^f^: 

M h 
3=1 

If the velocity field v n is continuous in fl, so that its pointwise values make sense, then we define 

Xj being the center of E^. Otherwise, assuming v n € (i 1 (fi)) , we take the approximate velocity 
field u\ to be defined by 

U !h = j V n{x)dx, 
E h 

which, owing to Lebesgue differentiation theorem, converges to v n in (L 1 (f2)) d as the grid is refined, 
i.e., when h —> + . Notice that in both cases u\ converges pointwise to v n for h —> + (up to 
possibly passing to subsequences and discarding Lebesgue negligible subsets in the second case). 
Besides this approximation property, we further require: 

Assumption 2. Let h, At > be such that 

At max \\v% h \\<h, (2.17) 

j=l,...,M h " 



where \\ ■ \\ is any norm in 



Equation (2.171 reminds of the classical CFL condition arising for the stability of numerical 
schemes designed to approximate hyperbolic conservation laws. In the present context, it imposes 
a limitation on the maximum displacement of any grid cell E^ produced by g 1 ^, so that when 
computing 

C d {{9t)- l {E*)) = C d {{9tT l {E>!) n E)) = C d {El n 

3=1 3 = 1 

we see that there is a fixed maximum number, say c, of non-empty intersections E^ n g^(E^) for 
i, j running from 1 to Mh, hence 

C d {{gl)-\E^)<cC d {E^) = ch d . 
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In addition, g™ L being a piecewise translation on fl, the measure g^jfcC is absolutely continuous 
with respect to C d , and its density r% £ satisfies, for a.e. x £ O, ^(x) > and: 

For R sufficiently small and for a.e. x £ ft the ball Br(x) is fully contained into one of the grid 
cells, thus in view of Eq. (2.17) we have 

, ,. c£ d (B R (x)) 
r?(x) < lim „ ./ ; // = c 
hK > ~ R^o C d (B R (xj) 

and finally 

^{{girHE)) = J rt{x) dx < cL d (E), V£e B(fi), (2.18) 

E 

which compares with the analogous property of 7„ stated by Assumption [T] 

Bearing all of this in mind, a numerical scheme is obtained by imposing that the P^'s satisfy 
Eq. (2.15) with the choice rj — Xe h i * = 1> ••• > Mh, under the action of the mapping g^: This 
gives 

r pZ +1 (x)dx= J pj;{x)dx, 

whence, using the specific representation of P£, P^ +1 , the following scheme is derived: 

Pi* = p E Plh^glT 1 ^) n E% i = 1, . . . , M h and n > 0. (2.19) 

j'=i 

It is customary for the sequel to define the measures {A£}„>o such that 



d\ n h :=PJidx. 



Notice that Eq. (|2.19[) does not correspond to the push forward of AJJ by g?, i.e., A^ +1 ^ 5/^#A^, 



indeed the measure 3^#A^ is in general not piecewise constant on the grid {Ej}j^, even for piece- 
wise constant A? and for a piecewise translation while A^ is forced to be so by construction. 
However, the scheme (2.19) enjoys some nice properties that we briefly list below: 

(i) It is positivity-preserving, in the sense that, given a discretization P^ of the initial density 
po such that p® h > for all j, it results p" h > for all j and all n > 0. This is consistent 
with the expected nonnegativity of the true densities p n 's. 

(ii) It is conservative, indeed 

M h M h M h 

ii^ +1 iii = ^E^ 1 = E^E £d ^)" 1 (^) n£; i) 

i=l 3=1 i=l 

M h M h 

= J2pl h C d ((g^- 1 (Q)nE^ = h d J2plk = Will, 

J=l 3=1 

hence 

1^11! = KHx, Vn>0, 

in accordance with the analogous property stated by Theorem [3] for the p n 's. 

(iii) It is boundedness-preserving, in the sense that if the initial datum po is discretized in such 
a way that P° £ L°°(SY) uniformly in h, i.e., there exists a constant Bo > 0, independent 
of h, such that 

Halloo = . max \p° j h \ <B , 
then P£ £ L°°(fl) uniformly in h as well, with 

||P£IU < c n B Q , Vn>0. 
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This follows by induction from Eqs. (2.19) and (2.18), whence we get 

I -Pi 



pn+l II 

\ r h Moo 



max \p^\<. < 4P£\U 

M h ' t=i,...,M h nr 



and is again consistent with both the boundedness of the true densities p n claimed by 
Theorem [3] and the related estimate on their L°°-norm. 



Given any piecewise constant function Ph € L 1 (0) n L°°(Cl) on the grid {^}j = V 



AW = X! Pi,hXE*( x )i Pjh > V j = 1, .. . 



and defined the measure dXh = Ph dx, let us denote by g^Xh the measure whose density with 
respect to C d is the function 

M h ( M h \ 



Notice that, using this notation, we can rewrite the scheme (2.19) compactly as A^ +1 = gjJjJAJJ. 
From the previous reasonings we know that, in general, (g%$Xh){E) ^ (g%#Xh){E) for E £ B(Q), 
however: 

Lemma 4. For grid elements {E^} k ^ it results 

(9MX h )(E%) = {gl#X h ){E h k ). 

Proof. (i) We preliminarily observe that g^ifcXh <C C d , indeed let E £ B(£l) be such that 
C d (E) = 0, then 



( 9 U*h){E) = X h ((gZ)- 1 (E)) 



Ph(x) dx 



< \\PHUL\igir\E)) < \\P h \\ oo cC*{E) = 0. 

The density g%#P h £ L 1 (O) D £°° (Q) of 9h# x h with respect to C d is given, for a.e. x £ Q, 
by 



(g%#X h )(B R (x)) 



= EPi.hX fl -(B*)(a5)- 



j'=i 



(ii) Given any grid element F^, 1 < k < Mh, we compute: 

{g n JX h ){E h k ) = / (gZ$P h )(x) dx = £ n 



and we have the thesis. 



□ 



From this result we can derive some stability properties of the scheme. 



12 



BENEDETTO PICCOLI AND ANDREA TOSIN 



Theorem 5 (One-step stability). Let p : B(Sl) — > R + be a measure supported in SI, with density 
p G L 1 (Sl) n L°°(Sl) with respect to C d . Assume moreover that j n : SI SI is a diffeomorphism. 
Then, for h > sufficiently small, there exist two constants A, B > such that 

M h 

E |(7n#M)(4) ~ (g^\ h )(E'j)\ <A\\p- nllx + Bh. 

Proof. (i) We begin by observing that 

|( 7 „#/i)(4 l ) - (g^\ h )(E';)\ < |( 7 „#/i)(4) - ( 7 n#A fc )(4)| 

+ |( 7 n#A h )(4)-(3^A /t )(4)|, 

and further that, owing to Lemma [4j 

= |( 7 „#/i)(4)-(7n#A h )(4)| 

Set 

IiJ ~ |(7n#M)(^) - (7n#A h )(^)| , I 2J := |(7»#A/0(4 ) - G##Afc)(4 )| . 

(ii) It is an easy consequence of Theorem [2] that "fnftp is absolutely continuous with respect 
to £ d . Let 7„#p G L^Q) n L°°(Q) be its density, then 

( ln #p)(x)dx= J p(x)dx = y'p( 7 ^ 1 (y))|det( 7 - 1 )'(y)|^, V£e6(!]), 

£ ln X (B) E 

where ( 7 ^ 1 )' is the Jacobian matrix of 7~ 1 . From the arbitrariness of E we get then 
(7n#p)0) = P('7r7 1 (^))l det (7^ 1 )'(^)|, for a.e. x G SI. 

(hi) We estimate now Since 7n is a diffeomorphism, there exist two constants A + , A_ > 
such that 

||( 7 „)'(x)|| < A + , ||( 7 - 1 )'( a; )|| < A-, for all x G fl, 

whence 

Xi,i < / |p(7^ 1 (^))-^(7^ 1 W)| |det(7- 1 )'(x)| da; 



<Al / |p(a;) - Ph{x)\ | det 7^ (a;) | dx < AlA+ J \p(x) - P h (x)\ dx. 
and finally 

for A = AiA^_. 

(iv) The estimate on I2 ,j requires a bit trickier reasoning as it involves the comparison between 
the push forward operated by the true motion mapping j n and the one operated by its 
approximation g^ . By a procedure analogous to that used in the proof of Lemma [4] to 
compute the density of g^^Xh, it is straightforward to see that the density of 7 n #A/ 1 can 
be also written as 



(y n #Ph)(x) = J2plhX lni E^)\ det ( 7 - 1 )'(x)|, 

3=1 
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whence 



dx 



E'J 



< 



l^llooE 



i = l 



A d _L d ( 7n (E?)nE%\gZ(E?)) 



| det (j- 1 )'^) - 1| dx + C d {g n h {E*) n E* \ 7n (£*)) . 
7 n (£f)ng£(£f)nE 3 h 

Recall the formula det (7^ 1 )'(a;) = (det 7 ^(y)) 1 for x = J n (y) an d denote I € M dxd the 
identity matrix. We have 7^(2/) = I + Atv' n (y), with furthermore ||^(y)|| < L all y € fi for 
a certain constant L > (due to the smoothness of v n ), whence |trw^(y)| < L d (tiv' n (y) 
being the trace of v' n (y)) all y £ il and finally 

det 7 ; (y) = 1 + A* tr v' n (y) + o(At) {At -> 0). 

Thus we conclude, for At — > 0, detj' n (y) ~ 1 + if At for a constant X e R. It follows 
(det7;(y)) _1 - 1 - KAt, therefore, for At sufficiently small, |(det7;(y)) _1 - 1| < \K\At, 
which allows to specialize the previous estimate of I2 ,j as 

M h 

l%i < \\PhUC* J2 [£ d (-Yn(E?)Agl(Et) n £*) 
i=l 

+ At£ d ( ln (E?)ngZ{E?)nE*)], 

where C* = max{Al, \K\, 1} and A denotes the symmetric difference of two sets. Sum- 
ming over j further yields 

M h M h 

5> 2J < ||P fc |UC* 2 [C d {ln{Et)Ag%{E*)) + At£ d ( ln (E?) n ■ 

It is now simple to obtain 

L d { ln {E*) n < £ d (gZ(Eh) = h d , 

whereas in order to estimate the Lebesgue measure of 7n (E^)Agl(E^) we rely on the 
following argument. Given x £ E^, the maximum distance between the points 7n (:c) and 
g^{x) is estimated as \"f n {x) — 9h( x )\ = At\v n (x) — v n (xj)\ < LAth, therefore, by taking 
x G dE 1 *, we see that the set j^E!?) \ g^(E^) is fully contained into the gap between 
the (hyper)cube g^(E^) with edge size h and the concentric (hyper)cube with edge size 
h + 2LAth. This entails, for At small (say At < 1): 

£ d ( 7 „(^) \ gUEh) <(h + 2LAth) d - h d = 0{Ath d ). 
Conversely, the set g^(E^ l )\j n (E l i l ) is fully contained into the gap between the (hyper)cube 
9h(Ei) w ith edge size h and the concentric (hyper)cube with edge size h—2LAth (assuming 
At < (2L) _1 for consistency). Consequently: 

£ d (.9/"(^) \ 7n(4)) < h d - (h - 2LAth) d = 0(Ath d ), 
whence finally £ d ('f n (Ej)Ag%(Ej )) = 0(Ath d ), which allows to complete the estimate of 
Ej-^j as (consider that M h h d = C d {n)) 

M h 

J2 X U ^ \\P h \\ooC**At, 
3=1 

C** being a new cumulative constant independent of cither h, At or n. Considering that 



condition (2.17) forces the order of magnitude of the time step At to be dictated by the 
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order of the grid size h, the previous relation implies finally / 2 j < Bh for a suitable 
constant B > 0. 
(v) Collecting the results so far obtained we ultimately have 

M h M h 

E |(7n#M)(4) " (fffc»A h )(4)| < +l2,j) <A\\p- P h \\i + Bh 

3=1 i=i 

and the claim of the theorem follows. □ 

Theo rem |5| applied to the measures {p n }n>o arid {A^} n >o gives the stability of the numerical 
scheme (2.19) with respect to the push forward (2.7) at one time step: 



E 



M„+i(4) - X% +1 (Ej)\ < A\\ Pn P^lli + 



This amounts to controlling the total variation of the error p n +i — A^ l+1 over the grid cells (left- 
hand side) by the ^-distance of the densities of the same measures at the previous time step and 
the size of the grid (right-hand size). By slightly modifying the argument used in this proof we 
can also obtain a stability over several time steps. 

Theorem 6 (Multistep stability) . Let 7„ be like in Theorem [5| and assume that v n is uniformly 
bounded in time over tt. Then, for h > sufficiently small, there exist a constant A n > such 
that 



max 



M4) - X n h (E'l)\ < . max U (£*) - A£(M)I + A n h d . 



j=l,...,M h ' J ' ■•-■>■> j=l,...,M h 

Proof. We proceed inductively on n and consider 

l/W^) - \Z +1 (E*)\ < |( 7 „#Mn)(^ 1 ) ~ (7n#A^)(4)| 

+ \( 7n #Xl)(E';)-(g^\ n h )(E^\ 



(i) The second term at the right-hand side is estimated like Iqj in the proof of Theorem 
[5| In particular, owing to condition ( |2.17[ ) along with the uniform boundedness of the 
velocity v n on f2, only a finite constant number of intersections j n (E^)Ag^(E^) n Ej, 
j n (E^) n g%(El l ) n Ej are nonempty when i runs from 1 to Mh, whence, for h sufficiently 
small, 

|(7n#A£)(l#) - (gl#X]\){E^)\ < \\PR\\ 00 Crh'» 1 , 

where C* is a constant independent of both h and n, and \\Pfi ||oo < c n Bo. 

(ii) The first term at the right-hand side coincides with Ii j in the proof of Theorem [5j but 
we need here a different estimate for it. Specifically, we write 

|( 7 „#M„)(^) - ( 7 „#A£)(4 l )| < | Mn (^) - \ n h {E^)\ 

\ Pn (x)-PZ(x)\dx, 

7^ 1 (Bj l )AB J ' 1 

then estimate 

\ Pn {x) - PR(x)\ dx < (C n \\p \\oo + c^Bo^-^E^AE*). 

The maximum distance between a point x € f2 and its image J n { x ) is 

\x- ln (x)\ < At\v n (x)\ < AtV, 



where V > 



comes for the uniform in time boundedness of the velocity v n on f2. 



Therefore, if At is chosen in such a way that AtV < h then condition ( 2.17 1 is fulfilled and 
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in addition, thinking now of a; € dj n the inverse image j n l (Ej) is fully contained 

into the (hyper)cube centered at Xj (the center of E 1 *) with edge size h + 2AtV, whence 

£ d (7« \ E'J) <(h+ 2AtV) d - h d = 0(h d ). 
On the other hand, we have 

\ % X {E))) < C d (E'l) = h d 

so that finally £ d (7~ 1 (i?j l )A£ , j l ) = 0(h d ) and for h small this term dominates over the 
previous one. 

(iii) Owing to the above estimates, we can find a constant a n > 0, which essentially depends 
on 11 through the powers of C and c, such that 

- X n h +1 {E*)\ < \p n (E$) - AS(£#)| + a n h d , 

hence taking the maximum over j = 1, . . . , and then summing telescopically over n 
we get the thesis with A n = Y^=o ai - ^ 

The quantities \p n (Ej) — A^(£^)|, j = 1, . . . , M^, can be interpreted as the localization error 
of the measure p n produced by the measure AJJ over the grid cells; in other words, they say how 
much the cell-by-cell distribution of the measure AJJ is a reliable approximation of the distribution 
of p n . We incidentally notice that the cell- by-cell distribution of AJJ is actually the most visible 
outcome when looking at the result of a numerical simulation. Theorem [5] guarantees stability, in 
a single time step, of the total variation of the localization error, but requires for this the stronger 
information on the Z^-distance between the exact and the approximate densities. Theorem [6] 
yields instead the uniform stability of the localization error over several time steps with respect 
to the grid size, starting from the initial error produced by the discretization of po. 

2.3. Nonlinear fluxes. In this subsection we briefly comment on how the previous theory applies 
when the velocity field v n depends explicitly on the measure p n , v n = v n [p n ], which gives rise 
to a nonlinear flux p n v n [p n ]. As a matter of fact, this structure turns out to be quite relevant 
for a wide range of applications, including pedestrian flow that we will address later, where the 
dynamics is strongly influenced by the current configuration of the system. 

Given the density po € L 1 ^), p > a.e. in ft, of p with respect to C d , let us focus on the 
existence of the densities {p n }n>o in case of motion mappings j n = 7n[yU n ] explicitly depending 
on the PnS. The first step consists in assuming that, for a certain n > 0, the measure p n is 
absolutely continuous with respect to C d , with density p n £ L 1 ^), p n > a.e. in ft. Next one 
checks whether, in view of this hypothesis, property ^ of Assumption [I] is actually satisfied. If 
this is the case, then Theorem [5] implies p n +i *C C d with density p n +i £ L 1 ^), p n +i > a.e. in 
f2, thus one concludes by induction that p n -C C d for each n > as desired. In practice, this is 
nothing but the same reasoning scheme applied in Subsect. |2. 1| to prove Theorem [2] with the only 
difference that property ^ of Assumption [l] on j n is not assumed to hold a priori, but is checked 
a posteriori as a consequence of the inductive hypothesis p n <C C d . Notice that uniqueness of 
the p n 's along with L 1 and possibly L°° estimates now follow straightforwardly, because the proof 
of Theorem [3] doe s not rely on any specific structure of the motion mapping 7„, except that it 



satisfies Eq. (2.13) 



As far as the spatial approximation of the jit n 's is concerned, one has to take into account that 
the approximate motion mapping g 1 ^ now depends on A^ because so does the approximate velocity 
field u"%. Basically, one tries to mimic the true velocity v n [p n ] by u£[AJJ], setting 

ulu [K] = v n [A£] (ay ) or [AJJ] - ^ f v n [A£] (x) dx 



and then reconstructing 



ul[\l]{x) = Y J ul h [K]xE^)- 
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By inspecting the proof of Theorems [5j [6] we see that a relevant part of the error estimate 
involving I2 relies on the Lipschitz estimate of the velocity v n within each grid cell Ej\ 

\ ln {x) - gl{x)\ <LAth, VxeE%. 

When the motion mappings 7„, g% depend on the measures /x n , AJJ, some considerations on the 
localization error produced by the scheme in a single time step might still be developed if the 
velocity v n is such that 

\vM{x) - v n [\ n h ]{x )\ = O (h+\\ Pn - PZWv^) , Vx G E). (2.20) 

We refrain from going into this issue, we simply point out that, as we will see more in detail in 
the application to pedestrian flow (cf. Sect. 3.2), an estimate of this form is more easily proved 
when 7„, g% depend on /i„, AJJ through p n , P£, respectively, in a non local way. 



3. Macroscopic modeling of pedestrian flow 

In this section we address the application of pedestrian flow using the time-evolving measure 
framework previously developed. Pedestrians in motion within a given walking area are an essen- 
tially discrete system, in which each person plays the role of an isolated agent moving in interaction 
with other surrounding agents, usually with the aim of reaching a particular destination. As such, 
this system should be conceptually described at the microscopic scale, i.e., by a system of ordinary 
differential equations tracing the evolution in time of the position of each single pedestrian. The 
possible coupling of these ODEs should reflect the interactions among pedestrians, namely the 
influence that each one of them has on the motion of the others. Some microscopic models of 
pedestrian flow are currently available in the literature, generally regarding pedestrians as rigid 
particles, whose motion is formally regulated by classical Newtonian laws of point mechanics: 

»i(t) = Fi{t), 

where Vi is the velocity of the i-th pedestrian, Fi an overall force acting on her/him, while the 
index i runs from 1 to the total number of pedestrian considered in the model. However, the force 
field Fi usually invokes nonclassical dynamic concepts, resorting mainly to the ideas of preferred 
direction of motion and of comfort/ discomfort due to the distance from or the proximity to other 
pedestrians. Therefore, the F/s have not to be understood strictly as mechanical actions exerted 
by pedestrians on each other. For instance, Helbing and coworkers [HHH] introduce the concept of 
social (or behavioral) force, which measures the internal motivation of the individuals in performing 
certain movements. Specifically, in their model pedestrians are regarded as points, and two main 
factors contribute to the definition of the social force acting on the i-th individual: 

(i) A relaxation toward a desired velocity, i.e., the velocity that the i-th individual should 
possess in order to reach her/his destination as comfortably as possible. 

(ii) A repulsive effect from neighboring pedestrians located too close, i.e., within the so-called 
private sphere of the i-th pedestrian, or from walls and borders found in the walking area. 

Essentially the same guidelines underlie the microscopic model of pedestrian flow by Maury and 
Venel [15]. The only difference is that pedestrians are now regarded as rigid disks, hence the 
repulsive effect among them is meant here as a geometrical constraint to avoid that disks step over 
one another. An alternative formulation, still relying on a microscopic description of the system, 
is instead used by Hoogendoorn and Bovy [111 112) , who propose a theory of pedestrian behavior 
based on the concepts of walking task and walking cost. Basically, they assume that pedestrians 
are feedback-oriented controllers, who plan their movements on the basis of some predictions they 
make on the behavior of the other individuals. Predictions are dictated by a sort of cooperative 
or non-cooperative game theory; in either case, they are affected by a limited in time and space 
predictive ability of the walkers. Each pedestrian behaves so as to minimize her/his individual 
estimated walking cost, which is expressed by a suitable functional depending on the predicted 
positions of the other people. 

A different approach to the description of the system uses instead partial differential equations 
and the theory of (possibly multidimensional) conservation laws. In this case, it is assumed that 
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pedestrians moving within a given walking area have a continuous distribution in space, so that 
it makes sense to introduce their density p — p(t, x) and to invoke some conservation principles, 
e.g. the conservation of mass and possibly also of linear momentum, in order to get an equation, 
or a system of equations, satisfied by p. Colombo and Rosini [5] propose a one-dimensional model 
built on a Cauchy problem for the nonlinear hyperbolic conservation law 

d t p + d x f(p)=0, xeR, t>0, 

which reminds of the Lighthill-Whitham-Richards (LWR) model of vehicular traffic (see Lighthill 
and Whitham [14], Richards [E]). The main difference is that the density of pedestrians exhibits 
two characteristic maximal values R, R* , with < R < R*, at both of which the flux / vanishes. 
In normal situations p ranges in the interval [0, R] 7 where the flux is nonnegative, either strictly 
concave or with at most one inflection point, and has precisely one local maximum point. When 
the density grows above its "standard" maximum value, i.e., for p £ (R, R*], it is assumed that, 
unlike vehicular traffic, pedestrians can still move but feel overcompressed, hence their flux is 
less effective than before and they enter a panic state. In this region, the function / features 
a trend similar to that described for p 6 [0, R], but with a local maximum value strictly less 
than the previous one. As illustrated in [SJ |B], this allows to define a concept of solution to the 
above conservation law in which non-classical shocks are admitted, i.e., shocks complying with 
the Rankine-Hugoniot condition but possibly violating entropy criteria. As a consequence, the 
classical maximum principle for nonlinear hyperbolic equations, stating that the solution p(t, x) 
remains confined within the same lower and upper bounds of the initial datum for all x G R and 
all t > 0, no longer holds true and the model is able to describe the transition of pedestrians 
to panic even starting from an initial density entirely bounded below the standard maximum R. 
The resulting fundamental diagram, i.e., the mapping p \— > f(p), agrees well with experimental 
observations reported by Helbing et al. in [8]. 

Bellomo and Dogbe [2] refer instead to a two-dimensional setting, in which the walking area 
is represented by a bounded domain O C R 2 with possible inlet and outlet regions along the 
boundary dil. The motion of pedestrians is described by a system of two partial differential 
equations invoking the conservation of mass and the balance of linear momentum: 

d t p + V • (pv) = 
d t v + (v ■ V)v — F[p, v], 

where F is a material model for the acceleration of the individuals depending in general in a 
functional way on the density p and the velocity v. The above equations are formally inspired by 
the classical fluid dynamics models of continuum mechanics, however the force (per unit mass) F 
contains non-classical contributions accounting for: 

(i) A relaxation toward a desired velocity, that makes pedestrians point in the direction of a 
certain target they want to reach. 

(ii) A local crowding estimate, based on the pointwise values of V/? possibly taken along the 
direction of the desired velocity, which might induce pedestrians to deviate from their 
preferred path in order to avoid areas of high density. 

(iii) A pressure-like term, possibly regarded as a material quantity as in the celebrated Aw- 
Rascle model of vehicular traffic pQ, which models the reaction of pedestrians to the 
presence of other individuals in the surrounding environment. 

Additional topics, like e.g. the existence of a limited visibility zone for each pedestrian when trying 
to evaluate the minimal crowding direction, are also discussed. In particular, special attention is 
paid to the characterization of the panic state and to the transition to it from regular conditions: 
The Authors of [5] suggest that pedestrians entering a panic state tend to follow chaotically other 
individuals, dropping any specific target, and therefore are mostly attracted toward areas of high 
density rather than seeking the less congested paths. 

As anticipated in the Introduction, the usually large amount of coupled ODEs to be handled 
simultaneously required by microscopic models is a drawback if not from the computational point 
of view, thanks to the increasing power of modern calculators, certainly for analytical purposes. 
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Among others, we like to mention here the chance to recover a global overview of the system 
from the knowledge of its microscopic state, possibly in connection with control and optimization 
issues. Macroscopic models are more suited to this, but those currently available in the literature 
have to face several other difficulties due to their intrinsic hyperbolic nature. First of all, the 
most natural settings for pedestrian flow problems are two-dimensional: One-dimensional models 
are essentially explorative, but they are unlikely to provide effective mathematical tools to deal 
with real applications. However, it is well known that the theory of multi-dimensional systems 
of nonlinear hyperbolic equations is much more complicated under both the analytical and the 
numerical point of view, therefore a sound mathematical mastery of such models may hardly be 
achieved. Secondly, the imposition of boundary conditions may be tricky in hyperbolic models, 
because on the one hand one is forcedly driven by the characteristic velocities in defining the 
inflow and the outflow portion of the boundary, while on the other hand it must be guaranteed 
that pedestrians do not enter or exit the domain from any point of the boundary other than the 
prescribed inlet and outlet regions. This issue gets even more complicated in presence of obstacles, 
which have to be understood as internal boundaries to the walking area. We notice, however, that 
the most interesting problems for applications generally do not concern pedestrian motion in free 
spaces, but precisely in areas scattered with obstacles (e.g., pillars, bottlenecks, narrow passages, 
see Helbing et al. [TO]); sometimes used to force the flow of crowds in specific directions. 

A time-evolving measures approach to the modeling of pedestrian flow may help overcome 
some of the difficulties just outlined. If one understands the mapping fi n as a mathematical tool 
to evaluate the degree of space occupancy by pedestrians, then it is perfectly natural to address 
the problem from a macroscopic, even Eulerian, point of view in spite of the intrinsic Lagrangian 
granularity of the system: Given a certain walking area Q C R 2 , the number /jL n (E) is simply a 
measure of the crowding of the subset E C fl at time n. In addition, as shown by the discussion 
of Sect. [2| there are basically no differences in the one- or two-dimensional theory, therefore 
one can immediately tackle realistic problems without the need for conceiving preliminary one- 
dimensional approximations of them. Finally, the imposition of at least the most common no-flux 
boundary condition, including hence the treatment of internal obstacles to the walking area, is 
quite straightforward, for in principle it simply requires to guarantee that supp /i„ C Q for all n. 
Notice that if O C £1 is, e.g., an obstacle, then it might be sufficient to impose ^ n (0) = 0, which 
translates the idea that pedestrians cannot occupy the area covered by O. Then any other subset 
E of £1 is automatically measured by taking the presence of the obstacle into account, indeed, 
assuming conveniently that /i„ is complete, it results \i n (E fl O) = 0, thus /J, n (E) = pi n {E\0). 



3.1. Nonlocal flux. The mathematical structure depicted by Eqs. (2.7), (2.8) requires as major 
modeling task to devise the velocity v n at the n-th time step as a function of x € H and possibly also 
of the current distribution ji n of pedestrians. Following some ideas proposed by Maury and Venel 
15 in the frame of microscopic models of crowd motion, we distinguish two main contributions 
to the overall velocity of pedestrians: 

V n [Hn](x) = V d (x) + V n [ll n ](x). 

The desired velocity Vd = Vd(x) is the velocity a pedestrian would have in the absence of other 
surrounding pedestrians. This component of the total velocity describes the preferred direction of 
motion toward specific targets, possibly taking into account the presence of intermediate obstacles 
within the domain Q. Therefore, it is not affected by the actual crowding of the environments, 
but it specifically depends on the geometry of the walking area (in this sense, it is a sort of field 
velocity). 

Assume the target of pedestrians is a certain location xq € f2 representing, e.g., an aggregation 
point inside the walking area or a door along its boundary (as illustrated in Fig. [I]). The simplest 
case is when fl is star-shaped with respect to xq, so that from any point x € £1 there exists a 
straight path to Xq. Then the desired velocity at x is readily defined as a vector pointing to xq, 
namely 

v d {x) =p d (x)^ — (3.1) 
\xo - x\ 
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FIGURE 1 . The desired velocity Vd toward the target xq and the interaction neigh- 
borhood Br(x) for a pedestrian located in x e fi. 



where pd : Q — > R + is a scalar nonnegative function representing the magnitude of Vd (in other 
words, the desired speed of pedestrians). 

When the domain f2 is scattered with obstacles, one has to take into account that some points 
x £ cannot be directly connected to xq by straight paths. For such points Eq. (3.1| must be 



duly modified, for instance by identifying several intermediate targets to be preliminarily reached, 
so as to bypass obstacles before pointing to xq (see Fig. [2]). The definition of these intermediate 
targets is clearly not unique, and might be related to some strategy of walking optimization, for 
instance the minimization of the total length of the path from a; to so- The interested reader is 
referred to Maury and Venel [TS] for further details. In this early application to pedestrian flow we 
refrain from going too much into this issue, deferring to a forthcoming work a deeper investigation 
of this modeling aspect. 

The interaction velocity v n — v> n [fj, n ](x) expresses the deviation of pedestrians from their pre- 
ferred path due to the presence of other surrounding pedestrians. In particular, we assume that 
interactions among pedestrians are effective only within a finite distance R > 0, defining a neigh- 
borhood of interaction Br(x) for each point x G fi: 



B R {x) = {yeil 



<R}. 



The specific shape of this neighborhood depends on the M 2 -metric used in the above definition. 
If, for instance, || • || is the Euclidean norm then Br(x) is the closed ball centered at x with radius 
R (see Fig. [IJ. Instead, by taking || • || to be the oo-norm the neighborhood of interaction results 
in a square centered at x with side 2R. 

Resorting to the idea of nonlocal interactions among agents proposed by Canuto et al. j3] for 
rendez-vous problems, we assume that pedestrians interact in a nonlocal way with all subjects 
comprised within their neighborhood of interaction. The main difference with respect to the 
above-cited rendez-vous is that in this case, at least under normal, i.e., no panic conditions, 
instead of aiming at aggregating they carefully try to bypass highly congested surrounding areas. 
We propose two different ways in which this behavior can be modeled, both based on the detection, 
by a weighted average, of a suitable point x* € Bn(x). In the first case, x* is the center of mass 
of pedestrians in the neighborhood Br{x): 



fi n {B R {x)) J 
n 
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» -'0 



Figure 2 . A point x € ft that cannot be directly linked to the target xq £ dfl due 
to the presence of an intermediate obstacle O C ft. In this case, an intermediate 
target x' , coinciding with a corner of O, is identified, whence the actual target is 
then accessible. The desired velocity at x points to x' , while the desired velocity 
at x' finally points to Xq. 



which can be thought of as an indicator of the area of highest crowd concentration within Br(x). 
Pedestrians trying to steer clear of this zone are driven away from x* , hence their interaction 
velocity can be given the form: 



V n [p n ]{x) = p v [n n }(x)(x - X*) 
( 1 
(J-n(B R (x)) 



Pu[^n]{x) — , 1 , rr / (x - v)Xb r (_x) (y) d(i n {y), (3.2) 



Pi/[/i n ] : ft — > R + being a nonnegative function related to the intensity of the interaction. The 
superposition of Vd{x) and v n [fi n ](x) may finally lead pedestrians to deviate from their preferred 
path in order to avoid local jams (see Fig. [3| left). 

In the second case, x* is instead a point of low crowding within the neighborhood of interaction, 
which can be computed as a sort of "inverse center of mass". In practice, it should average out 
the neighboring locations y € Br(x), emphasizing those with low pedestrian density p n {y)- 

yXB R (x){y){u o p n )(y) dy 



XB R (x) 

(y)(uj o Pn )(y) dy 

where to : M+ — s- M. + is a nonnegative nonincreasing weight function. Once this point has been 
defined, pedestrians locally deviate toward it, hence we put 

V n [Pn]{x) = p v \pL n ]{x)(x* - x) 

(y - x)xB R (x){y)(u ° Pn){y) dy 

Pu[Pn]{x)- 7 . (3.3) 

xb r {x){v){u ° Pn){y) dy 
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FIGURE 3. Construction of the velocity v n [p n ](x) when x* is the center of mass 
of pedestrians (left) or a point of low crowding (right) in Br(x). 



Notice that, unlike Eq. (3.2), in Eq. (3.3) the dependence of v n on p n is achieved through p n , 
which explicitly requires that a density for the measure \i n be defined with respect to L 2 for all 
n > 0. Again, the superposition of v ( i(x) and Vnlprijix) make pedestrians modify their direction 
of motion looking for uncrowded surrounding areas (see Fig. [3] right). 



In the sequel we will concentrate on the form of v n [^n\ provided by Eq. (3.3), the case of 
Eq. (3.2) requiring, at least for /i„ <C C 2 , essentially the same calculations up to the use of the 



boundedness of the density p n claimed by Theorem [3j In more detail, the case that we will be 
able to analyze theoretically is when 

p d (x) = a\x - x\, p v [fj, n ](x) = p J XB R (x)(y)(u°Pn)(y)dy, 



for a, P > 0, which gives rise to 



v n \p n ](x) = a(x Q -x)+P I (y - x)xb r {x){v){u ° Pn)(y) dy. 



(3.4) 



As a consequence, the motion mapping j n depends on the measure p n , therefore the theory 
presented in Sect. [2] applies to the present model provided the discussion made in Subsect. |2.3| 
holds for the velocity w„[/x ra ] resulting from Eq. (3.4). 

Let us prove that 7„ complies with property ([2j) of Assumption [l] by showing that v n [fi„] is 
Lipschitz continuous in Q. This is self-evident for the desired velocity Vd{x) — a(xo — x), while, 
as far as the interaction velocity f„[/i n ] is concerned, we compute, for Xi, x% € fi, 



\v n \pn}{x 2 ) ~ y n [^n]{xi)\ = P 



y (xB R ( Xl )(y) - XB H (x 2 )(y)) (^ o Pn )( y ) dy 



+ J (x2Xb b (x 2 ){v) - %iXb r ( x1 ){v)) {u p n )(y)dy 

h 



< /?w(0) 



Co / | XB R (x 1 ) (y) XB r (x 2 ){y)\dy 
n 



\X2XBr (y) (y)\ dy 



where Cq > is a constant related to the boundedness of f2, such that < Cq for all x £ f2, 
and where we have used w(s) < w(0) for all s > 0, together with the inductive hypothesis p n > 
a.e. in f2, the weight ui being nonincreasing on [0, +oo). Notice that for any two sets A, B and 
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any two real numbers a, b it results 



\aXA(x) - bxB(x)\ 



if x G A \ B 
if x e B \ A 
b\ ifieinfi 

otherwise, 



thus in particular, when a = b = 1, \xa{x) — xb(x)\ — Xaab{x), where AAB = (A\ B) U (B \ A) 
is the symmetric difference of the sets A and B. Therefore we can carry the previous computation 
on as 

\v n [Pn]{x 2 ) - v n \nn]{x-i)\ < f3u J {0)[C n C 2 (B R (x 1 )AB R (x 2 )) + \x 2 \C 2 (B R (x 2 ) \ B R { Xl )) 

+ \ Xl \C\B R {xx) \ B R {x 2 )) + \x 2 - Xl \C 2 (B R ( Xl ) n B R {x 2 ))]. 

Assuming that the neighborhoods of interaction are closed balls with respect to the Euclidean 
metric in M 2 and invoking the invariance of the Lebesgue measure under translations, it is basically 
a plane geometry task to see that 



C 2 {B R {x l )AB R {x 2 )) is 

thus globally 



= 2nR 2 if \x 2 -xt\ > 2R 

< AR\x 2 - xi\ if |af 2 ~x\\ < 2R, 



C 2 {B R {x 1 )AB R {x 2 )) < 4R\x 2 - Xl \. (3.5) 
From this we deduce immediately 

C 2 (B R ( Xl ) \ B R {x 2 )) = C 2 (B R (x 2 ) \ B R {x x )) = ~£ 2 (B R ( Xl )AB R (x 2 )) < 2R\x 2 - x x \, 

whence finally, using further C 2 (B R (xi) n B R (x 2 j) < ttR 2 , 

\v n [Hn]{x 2 ) - v n [Hn]{xi)\ < /3oj(0)R(8Cn + irR)\x 2 - xi\ 

which yields the Lipschitz continuity of the interaction velocity of pedestrians. Analogous calcu- 
lations can be repeated for square-shaped neighborhoods of interaction. 

Owing to the above results, we conclude that v n [[i n ] is Lipschitz continuous in Q as it is the 
algebraic sum of two Lipschitz continuous functions. The reasonings of Subsect. |2.1| entail then 
that 7„ complies with Eq. (2.13), hence in view of Theorems [2j [3] we have proved: 

Proposition 7. Let fi <C C 2 with density p G n L°°(fl), po > a.e. in ft, and let the 

velocity field (3.4) be given. Then: 

(i) fi n <C C 2 for all n > 0. 

(ii) There exists a unique sequence {p n }n>o C L 1 ^) n L°°(f2), p n > a.e. in Q, such that 
dp n = p n dx each n> 0, solving Eq. (2.7) (or, equivalently, Eq. (2.15)) with po as initial 
datum. 

We can therefore speak of density of pedestrians (in the sense of Radon-Nikodym theorem), 
and recover the measure p n as 

p n {E) = / Pn {x)dx, VEeB(n) 



each n > 0. Remember that the /i„-measure of the subsets of f2, more than the pointwise values of 
the density p n , is in principle the actual meaningful macroscopic information describing the space 
occupancy by the crowd. 
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Figure 4. The grid cell E- with side h and its center Xi 



3.2. Numerical approximation. Spatial approximation of Problem (2.16) can be carried out 
by fixing preliminarily the walking area to be, for instance, the square f2 = [0, 1] x [0, 1] C K 2 . Let 
us partition f2 by a uniform orthogonal grid, made of square elements of side h, that we denote 
by {Eij}f*j =1 (from now on we drop the index h for notational convenience). In this specific 
two-dimensional application, it is useful to identify each element of the grid by a couple of indices 
i, j, which locate the center iy of the cell as 

'2t-l, 2j-l. 

-h, 



i, J = 1, 



M, 



2 ' 2 

see Fig. |4j and to deal consequently with a double-indexed numerical density P r ' 



A I 



P"(x) = PijXE^x), n>0. 



In particular, the numerical scheme (2.19) rewrites as 



1 



J2 pLr 2 (^-n ff "[A»](^ m )), 



M, n > 0, 



(3.6) 



/, m— 1 



where we have used the invariance of C 2 under the (piecewise) translation # n [A n ] to get 

£ 2 (( ff n [A B ])- 1 (B tf ) n E lm ) = C 2 (E tl ng n [X n ](E lm )). 

The dependence of the approximate motion mapping g n on the measure dX n — P n dx amounts 
naturally to a dependence of the measures C 2 (Eij D g n [\ n }(E[ m )) on the coefficients {py}fj = i 
of P n . Notice however that the form (3.6) of the scheme is more convenient for implementation 



purposes than the generic form (2.19), as it does not require to determine the inverse images of 
the grid cells with respect to the mapping g n . Rather, we incidentally notice that their images 
are straightforwardly obtained by translation as g n [X n ](Ei m ) = E[ m + u]^JA™]Ai. 



According to the discussion of Subsect. 2.3 g n complies with property (??) of Assumption [2] 
provided at each time step the CFL condition (2.17) is satisfied. In more detail, for every fixed 



pair of indices (i, j) the intersection D g n [X n ](Ei m ) is nonempty for at most nine adjacent pairs 
of indices (I, m), namely 

(I, m) = (t, j), (i ± 1, j), (i, j ± 1), (i ± 1, j ± 1), 
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C 2 {Eij nj„(%j+i)) UxAt 




FIGURE 5. The intersection between Eij and g n (Ei + i j+i), which is nonempty 
only if the horizontal and vertical components U\, U2 of the velocity J+1 are 
both negative. In this case, the Lebesgue measure (area) of the intersection is 
C 2 (E l3 n g n (E i+ltj+1 )) = U^U^At 2 . 



so that, denoting by U\, XJ% the horizontal and vertical component, respectively, of uj^J.Yjj, the 
coefficients C?{E%j n g n [X n ](Ei m )) of the above scheme can be duly computed as (cf. Fig. p3J) 



C 2 {E l} n g n [\ n ](E lm )) = [U+AtS^ + (h - I^A*)^ + C/f At<5 M+1 ] 

x p^AtSmj-i + (h- \U 2 \At)S mj + U^AtS mJ+1 }. 



In this formula, S rs is the Kronecker delta: 



S rfi 



1 if r = s 
if r ^ s, 



whereas 



U + = max (U, 0), U~ = max (-17, 0) 



are the positive and negative part, respectively, of U . 

With respect to the discussion proposed in Subsect. |2.3[ let us compute, from the definition of 



v n given by Eq. (3.4) 



v n [Hn]{x) - v n \\ n ](C) =a(£-x)+0 / (y - x)xb Pn){y) dy 



P J (y-OXB Rii )(y)^oP n )(y)dy 
h 
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whence, adding and subtracting repeatedly several terms, 



\v n [[i n }(x)-v n [\ n }(0\ <a|£-x| + /3 



(y - x){xb r (x){v) - XB R {i){y)){u ° p n )(y) dy 



+ / (y - x)xB R (s)(y)(u ° p n - w ° P n ){y) dy 



n 



J XB R{i) (y)^°P n ){y)dy 

a 

< a\£ -x\ + 2(3CnH0)£ 2 (B R (x)AB R (O) + L u \\p n - P"^} 



where we have assumed that the weight function lo is Lipschitz continuous with Lipschitz constant 



> 0. Using Eq. (3.5 1 we further obtain 

< [a + 80C a u(O) + /3£ 2 (ftM0)]|£ -x\ + 2{3C Q L u \\p n - P n \\, 
=:C%-x\+C"\\p n -P n \\ l , 



which is precisely an estimate of the form (2.20) 



4. Conclusions and research perspectives 

In this paper we have addressed the modeling of pedestrian flow problems from a macroscopic 
point of view, resorting to classical measure theoretical methods and specifically to the concept of 
discrete-time evolving measures. 

The main idea consists in describing the space occupancy by pedestrians through a positive 
measure p n , parameterized by a discrete time index n £ N, which evolves according to the recursive 
push forward equation p, n +\ = 'YnifcPn- Here, "f n (x) = x+v n (x)At is the so-called motion mapping, 
expressing the space displacement during a time step At under the effect of a velocity field v n . 
Given a walking area O C K 2 , the number p, n (E) represents, for every measurable E C 17, a 
measure, in macroscopic averaged terms, of the crowding of the area E or, in other words, an 
estimate of the amount of people contained in E. Hence, the point of view on the system finally 
provided by the sequence of measures {p n }n>o is essentially Eulerian, in spite of the Lagrangian 
description of the motion yielded by the mappings {7n} n >o- 

The first part of the work (cf. Sect. [2]) has been devoted to a theoretical study, for generic 
d-dimensional systems, of the mathematical structures just outlined. In particular, we have first 
shown that the discrete-time push forward can be formally derived from an explicit time discretiza- 
tion of the classical mass conservation equation of continuum mechanics (cf. Proposition [lj. Then 
we have addressed specifically the question of the existence and uniqueness of a density p n for the 
measure p n with respect to the Lebesgue measure C d (cf. Theorems [2] [3J , so as to be able to speak 
of density of pedestrians at least in the sense of Radon-Nikodym theorem. Notice that if, on the 
one hand, measuring the space occupancy by p, n makes perfectly sense in a macroscopic frame, 
claiming as a modeling assumption that pedestrians can be physically described by a continuously 
distributed density may be, on the other hand, more hardly acceptable and more difficult to jus- 
tify. This issue is actually common to many other nonclassical systems, for which the concepts 
of mass and mass density are not as straightforward as in standard continuum mechanics. The 
theory we have developed in this paper provides, as a by-product, a possible methodology that, 
starting from reasonably meaningful macroscopic concepts, may help in gaining, at least in the 
abstract, a pointwise density for the system under consideration. 

The existence of a density is useful for the numerical treatment of the equations of the model, 
indeed approximating p n as a mapping on R d is undoubtedly more practical than approximating 
fi n as a set mapping defined on some a- algebra of M. d . In the paper we have proposed a numerical 
scheme to discretize, at each time step, the push forward of the p^s over a suitable space grid in 
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FIGURE 6. Desired velocity in presence of two obstacles Oi, C>2- The field is 
generated by a potential attracting pedestrians toward the target Xq. 

the domain O. In addition, we have studied the error produced by the numerical solution both in 
a single step of push forward and on the true density p n , on the basis of the level of refinement of 
the space grid compared to the discrete time step At. The scheme has proved to be robust and 
efficient in providing potentially accurate approximations of the true solution (cf. Theorems [5j 
[6]), as well as in preserving some interesting theoretical features of the latter, e.g., nonnegativity, 
integrability, and local in time boundedness. 

As mentioned above, the existence of a density is often also helpful in order to look conceptually 
at the system from a point of view close to that of classical continuum mechanics, even when, like 
in the case of pedestrian flow, the dynamics is mainly dictated by nonclassical evolution rules. 
In the second part of the work (cf. Sect. [3]) we have been concerned with the application of the 
time-evolving measures theory to the modeling of crowd motion. In particular, we have proposed 
a form of the velocity of pedestrians inspired by two main guidelines, somehow common to several 
microscopic models of pedestrian flow already available in the specialized literature: (i) The will of 
pedestrians to reach specific targets within the walking area, bypassing at the same time obstacles 
possibly present along their paths (desired velocity); (ii) The tendency of pedestrians, at least in 
normal, i.e. no panic, conditions, to look for uncrowded areas in their neighborhoods, in order 
not to stay too close to one another (interaction velocity). Notice that the simultaneous presence 
of Eulerian and Lagrangian aspects in the description of the system allows to duly incorporate 
microscopic viewpoints, namely the motion mapping j n and consequently the construction of the 
velocity v ni within a global macroscopic perspective on the system provided by the measures [i n . 

At present, our treatment of the desired velocity is mainly illustrative. In particular, we are 
assuming that pedestrians aim at one specific target, for instance an aggregation point or a door 
along the wall, that they can reach either directly or by possibly avoiding, in a simple manner, 
a single intermediate obstacle within the walking area (cf. Fig. [2]). However, it is plain that the 
most interesting problems concern the flow of people in areas scattered with obstacles, like e.g., 
pillars, bottlenecks, narrow passages, therefore we are currently scheduling a forthcoming work 
in which we will further investigate this modeling aspect, relying on the essentially geometrical 
nature of the desired velocity field. As an anticipation, we observe that it is for instance possible 
to generate smooth fields from a potential u : fl — > M, which attracts pedestrians toward a target 
xq G dtt, by solving the elliptic problem 
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where g vanishes on all external and internal boundaries of £1 but the one containing x , and then 
setting Vd(x) — Vti(x) (see Fig. |6|. 

Concerning the interaction velocity, the main characteristic is that it accounts for nonlocal 
effects within a certain finite-radius neighborhood around pedestrians. In practice, each individual 
is assumed to have her/his walking program influenced, in average, by all people comprised in some 
neighborhood of interaction, in such a way that she/he can possibly agree to deviate locally from 
her/his preferred path in order to bypass crowded areas, looking preferentially for free space. In 
the present work we have shown that a nonlocal average structure of the interaction velocity fits 
well in the proposed time-evolving measures theoretical framework (cf. Proposition [7]) . From the 
modeling side, this proposal deserves further investigation and development in order to provide a 
simple but accurate description of the interactions among pedestrians, which will be in turn the 
object of the above-announced subsequent work. 
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